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Effective constituent quark models d la Nambu Jona-Lasinio have recently been 
used to a large extent to describe the low-energy behaviour of QCD and the hadron 
spectrum (see [p]] for a recent review). Phenomenologically these models are found to 
be very successful in reproducing the experimental values of the low-energy coupling 
constants to 0{p i ) in Chiral Perturbation Theory (xPT) [@. Many of the couplings 
between resonances and pseudoscalar mesons have also been computed @, |3| and 
compare well with experiment, as well as the vector and axial- vector masses [0, |J. 

When we turn our attention to the scalar sector the situation is less satisfactory. 
Within the Extended Nambu Jona-Lasinio (ENJL) model the scalar two-point func- 
tion shows up a pole at a mass M$ = 2Mq Q, where Mq is the constituent quark 
mass. With a typical value of Mq = 250 -v- 350 MeV (this range of values is taken 
from the different fits of 0]) one has Ms = 500 4- 700 MeV. Experimentally there 
is no evidence for so low a scalar resonance. Until 1976 the Particle Data Book 
quoted a scalar resonance (called e) between 600 and 800 MeV with a width of more 
than 600 MeV. Because of this very large width it has then been discarded from 
the particle tables. Unclear experimental signal of a narrow scalar state around 750 
MeV is reported in ||, while the first clear scalar resonances are the ao(983) and 
the isosinglet /o(975) states. The proper qq assignment of the a (983) remains a 
problem. Its interpretation as ordinary qq system || is not excluded; other possible 
interpretations are as a qqqq state or a KK bound state 0. For the / (975) 
the analysis of ref. || supports the prediction of the unitarized quark model which 
interprets /o(975) as a qq resonance with a large admixture of KK virtual state, 
while ref. [|L(J almost excludes its interpretation as a KK molecule and favours 
a conventional Breit-Wigner-like structure. In both cases the interpretation of the 
/o(975) as an ordinary qq state is almost excluded. Recently, a fit of the available 
data has been performed [11] indicating that the KK component is large for both 
the a (983) and / (975) states. 

In the large- N c limit of the ENJL model the isoscalar resonance situated at 2Mq 
is predicted to have a very large width of more than 600 MeV [12). Assuming that 
such a resonance exists, qualitatively at least, there would be no surprise if l/N c - 
corrections in relation to the large isoscalar width created a non- degeneracy between 
the isoscalar (the hypothetical e) and the scalar octet (which could be the do). 

This is the physical intuition we would like to precise. The framework of the 
ENJL model offers a starting point to study the l/iV c -corrections. At this point it is 
worth mentioning that the NJL-like models contain most of the effective low-energy 
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models discussed in the literature as particular cases. Therefore such an approach 
lies in a quite general context. 

The effective ENJL Lagrangian is given by f2|: 

Cenjl = Cqcd + ~^r [(gg) 2 - (<?75?) 2 ] - jif KqinQ) 2 ~ (Ql^q) 2 ]- (1) 

The problem of the connection between QCD and this Lagrangian has been ad- 
dressed in [0]. The Lagrangian (|l|) is the first term of a double expansion in 1/N C 
and in 1/A^, A x ~ 1 GeV being the ultraviolet cutoff of the low-energy effective 
theory. 

The effective action derived in [0 from the Lagrangian (|I]) is subject to two 
types of 1/N C corrections. The first type will not be studied here. It concerns the 
1/N C expansion mentioned above: at next-to- leading order in 1/N C , there are other 
dimension six fermionic operators compatible with the symmetries of the original 
QCD Lagrangian. The second type appears when one considers loops made of chains 
of quark bubbles, i.e. loops of mesons. In this article we shall be interested in the 
latter. 

In section |I] we explain the relation between the quark-bubbles resummation in the 
non-bosonized ENJL model and the bosonized version. Then, in section |2| we present 
the effective Lagrangian we shall use in order to perform the loop calculation. The 
coupling constants of this Lagrangian are derived within the ENJL model. Section 
^] deals with the analysis of the l/A^-corrections. We show the corrections to the 
gap-equation, the mass-splitting between the scalar singlet and non-singlet, and the 
overall shift of the singlet. Finally section |4] is devoted to conclusions and remarks. 

1 Bosonized versus non-bosonized version of the 
ENJL model. 

Let us consider the isoscalar two-point function Il(g 2 ) —ij d 4 x e igx (0\TS(x)S(0)\0) 
where S(x) = —^q(x)q(x). This Green function has been computed in the large- N c 
limit within the non-bosonized version of the ENJL model using proper-time regu- 
larization with an ultra-violet cut-off A x . The result 0] is obtained by resumming 
the geometrical series 
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CXJ 

n( g2 )( ^oo) = n(g2) J2 (g s nm) n , (2) 



n=0 



where 

n(Q 2 ) = - - (Q 2 + (2M Q ) 2 )Z S (Q 2 ) (3) 
9s 

is the bare fermion loop diagram in the mean-field approximation, Zs{Q 2 ) is the 
wave function renormalization constant 

2 I ]\/T2 



, „ N N r r 1 ( a(l-a)Q 2 + MA\ 

z ^=t^ 2 1 dar (°-^4f — -) (4) 



with r(0,e) = / e °° dz\e~ z and g s = 4n 2 G s /N c A 2 x . One obtains 

(g) -i-^n ( g 2) - ^ Q 2 + (2M Q y)- (5) 

The next-to- leading contribution in 1/N C is built by dressing the leading- N c scalar 
two-point function with one loop of linear chains of constituent quark bubbles. We 
shall concentrate here on the example of the self-energy diagram, as shown in Fig.l, 
involving only the isoscalar coupling g$. 

These multiloop diagrams contain overlapping divergences. The regularization 
procedure we use is as follows: first we regularize the momenta running in the 
quark-bubbles and secondly we regularize the remaining momentum which runs in 
the internal loop made of quark-bubbles. This allows us to use the expressions found 
in (U for the two-point functions, like Eq. fl5|). The divergences occurring in the loop 
calculation have to be interpreted as physical divergences since the ENJL model is 
an effective non-renormalizable theory in which all the momenta are cutoff at the 
ultra-violet scale A x . 

Notice also that expressions like (^) should not be used beyond the convergence 
radius of the geometrical series, given in this case by q 2 = —Q 2 = (2Mq) 2 . Since 
2Mq is of the order of A x it is legitimate to make the integration over the internal 
loop, i.e. the bosonic loop, up to such energies. 

In order to compute the self-energy diagram of Fig. 1 we need the three-point 
function 
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f (q, k) = J d A x d A y e iqx e iqy (0\TS(x)S(y)S(0)\0). 



(6) 



It is found to be 



f(q,k) 



16tt 2 



4M C 



r o 



M, 



Q 



"ST' A: 



+ 0(^Y,(,-^). (7) 



We now proceed to the resummation ggU. = Yzffn + lZ^g Pa^ i-q<?n 



n(g 2 ) = — i 



[ZsmgsT 1 



with 



ffs V g 2 + (2M ) 2 - s(g 2 )z 5 (g 2 )-\ 



■ and get 

(8) 



l-^n(F) l-0 5 n((g-A;) 2 )- 



(9) 



This expression has to be compared with the one obtained within the bosonized 
approach, where we consider the three-scalars vertex at lowest order in the derivative 
expansion with coupling A: 



n(g 2 ) 



i 



(Zstygs)- 1 



where 



g s V g 2 + (2M Q )2-A 2 /(g 2 



d 4 k 



(10) 



(2tt) 4 F-(2M q ) 2 (g - k) 2 - (2M Q y 



Considering now Eqs. (||) and (U) and substituting T(g, fc) with T(0, 0) and Zg(g 2 ) 
with Z s {0), Eq. (|) reduces to Eq. (0) with 



T(0,0) 

z ff (°) 8/a ' 



(12) 
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This is nothing but the expression of the three-scalars vertex coupling A3/3! which 
is obtained in the next section (Eq. [15]), within the effective action approach []. 
Numerically the error made in approximating T(q, k) with T(q, 0) is of the order 
of thirty percent when k 2 is below A 2 . As a consequence of this approximation we 
shall not identify the cut-off used in the regularization of the internal loop, i.e. the 
bosonic loop, with A x , but rather keep it as an adjustable parameter whose order 
of magnitude is A ~ 2Mq ~ A x . 

In what follows we shall use the bosonized ENJL model where, for each vertex, 
only the lowest order in the derivative expansion will be kept. This approximation 
allows us to simplify the calculations and preserves chiral invariance. 



2 The effective Lagrangian 

We restrict ourselves to the SU(2) l x SU (2) # case. The general form for the meson 
matrices, singlets or triplets under SU(2)v, reads 

M = E^ M (^ (a) + ^ M ol, (13) 

where r^ a \ a = 1, ..3 are the Pauli matrices with Tr(r a T b ) = 25 ab and Mq is the sin- 
glet component. We use for the pseudoscalar field the exponential parametrization 
*(*) = e*p(-^£). 

In the chiral limit (m u = = 0), the effective chiral Lagrangian including scalar 
and pseudoscalar mesons to 0(p 2 ) is given by: 



£ S > P = I < > +1 < d » Sd " S > ~\ M S < S2 > +^nt 

££f = < s 3 > < > + Cd < su" > < s%e > +4 2) < s^se >, 

(14) 

where S is the scalar meson field and £ M = ^{^(^ — ir^)^ — £(<9 M — il^)^} is the 
axial current constructed with the pseudoscalar field, where r M and are the external 
right-handed and left-handed sources. The couplings among mesons are obtained 
by integrating out the constituent quark fields in the bosonized ENJL model using 
the heat kernel expansion with proper time regularization as explained in ref. 0]. 

^notice however that Z$ differs slightly from the expression given by the heat-kernel calculation; 
see [[§ about this question. 
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They are functions of the ultra-violet cut-off A x , the constituent quark mass Mq, 
the axial-pseudoscalar mixing parameter qa and the number of colours N c . We find: 



A3 
3! 

A4 
4! 

Cd 



r (i) 



X, | M, 



167T 2 

iV c l' 
16tt 2 Zf 



3/2 



r(o,e)--r(i,e) 



Nr. 



r(o,e) 

2^1 



4r(l,e) + -r(2,e) 



167r 2 ~'~'VZ s 
,2 



1 AT C ^ 

2 16vr 2 Z 5 

1 AT C g\ 

2 16tt 2 



r(o,e)-r(i,e) 



20, 



r(0,e)- y r(l,e) + -r(2,e) 

r(o, e )-Hr(i, e ) + lr(2,e) 



(15) 



The incomplete gamma functions T(n — 2, e), with e = Mq/A 2 , are defined as 

1 



T(n-2,e) 



z 



(16) 



r(— 2,e) contains a quartic divergence, T(— l,e) a quadratic one, T(0,e) is logarith- 
mically divergent, while T(n, e) with n > are finite. 

The Interaction Lagrangian of the scalar mesons with vectors and axial-vectors 
at the first two leading chiral orders, i.e. 0(p°) and 0(p 2 ), isfl: 



= ~c A < SA^A" > +c { l ] < SA^SA" > +cj } < S 2 A^ > 
+c A p < Situ, A"} > +4 ] < SV^SV^ > +4 2) < S%V» > 
+d v < SV^V^ > +d ( y ] < SV^SV^ > +d v V) < S 2 V^V» V > 
+d A < SA^A^ > +d { l ] < SA^SA» V > +d { l ] < S 2 A^ U A^ > . (17) 

Notice also the presence at 0(p°) of the mixed term scalar-pseudoscalar- axial with 
coupling cap- Zy and Za are the wave function renormalization constants of the 
vector and axial-vector fields: 

/V l 

Zv = ZA = Y^3 me) (18) 

The couplings are found to be: 

II The Lagrangian involving vectors starts at 0(p 2 ). For consistency we need to include 0(p°) 
as well as 0{p 2 ) terms in the Lagrangian involving axial-vectors. 
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ca 



C A 



(2) 
C A 



cap 
c v 

'I A 



4 } 



(i) 



d 



dP 



d 



(2) 



1 N, 



Mr 



2 16tt 2 ^fZTsZv 



4T(0,e)-4T(l,e) 



1 JV„ 



2 167T 2 

1 JV„ 



2 167T 2 



r(o,6)-^r(i,e) + ^r(2, e ) 
r(o,6)-^r(i,e) + ^r(2,6) 



1 N c g A M Q 

2 16tt 2 \[Z~ s jz v 



-4r(o, 



J2) 



1 JV„ 



y " 2 16vr 2 Z S Z V 

I N r 1 1 

-ri,e 



6 167T 2 Mg v^^V 

iV c l ' 1 



r(o,c) + -r(i,e) 



16tt 2 Mq v^^y 
1 iV c 1 1 
_ 2 16tt 2 Z^y-M^ 



"ir(i )C )-ir(2,e) 
|r(M) + lr(2,e3 



1 JVc 1 l_ 

'216tt 2 Z S Z V M% 

1 JVc 1 1_ 

'216tt 2 Z S Z V M 2 Q 
1 JVc 1 ]_ 

'2167T 2 Z 5 Zy M 2 



ir(i, f) + lr ( 2, £ ; 



1 



--r(i, e ) + -r(2, e ) 

1 T\ 
- 5 r(l,e) + -r(2,e) 



(19) 



3 Analysis of the 1/AT C corrections 



In the Appendix we have listed the results for all the loop contributions. We work 
in the chiral limit (m n = 0) and in the zero momentum limit (q 2 = 0). The masses 
of the vector and axial-vector mesons have been fixed to their experimental value: 
My ~ 0.770 GeV [13] and Ma — 1.2 GeV (we disregard the error on the axial mass). 

We have calculated the 1/N C corrections in two cases: 1) assuming that the scalar 
particle is a singlet (which, as we said in the introduction, cannot be identified with 
the / (975) resonance) 2) assuming that the quark content of the scalar particle is 
the same as that of the p(770) vector meson. This could be the case of the physical 
ao(983) scalar resonance. Obviously our 577(2) calculation has to be interpreted as 
a first indicative approximation of the fully realistic SU(3) calculation. 

The diagrams which contribute to the scalar two-point function at next-to-leading 
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order in 1/N C are shown in Fig. 2. They are the self-energy (a), the tadpole (b) and 
the 'top' (c) diagrams. 

Particles running in the loop are the singlet and the £77(2) triplet states for each 
quantum number: Scalar, Pseudoscalar, Vector, and Axial- Vector. We disregard in 
the present analysis the splitting between the singlet and the triplet pseudoscalar 
states; this is due to the £7(1) axial anomaly, which appears in the effective La- 
grangian at next-to- leading order in the 1/N C expansion. 

The correction to the mass depends quartically on the UV cut-off A of the bosonic 
loop. It is proportional to 



A V 



N c \2M Q J 

Hence reasonable values of A cannot exceed 2Mq. 



(20) 



Our conclusions are threefold: 

1) The top diagram can be alternatively included in the corrections to the scalar 
propagator or in the gap-equation so that they translate into corrections to the 
leading- A^ c value of Mq. These corrections do not create any mass splitting between 
singlet and non-singlet states. In addition they are q 2 independent. The main 
result concerning these diagrams is the compensation between negative contributions 
(pseudoscalar and scalar) and positive contributions (vector and axial) so that the 
corrections to the gap-equation are small. 



In Ref. [[L4[ an attempt to compute the l/A^-corrections to the gap equation of 
the NJL model has been performed in the Gy = case. However, the formalism 
they have adopted implies non- derivative couplings of pseudoscalar mesons to other 
particles. This makes the comparison with our approach problematic. In particular, 
the relative size and sign of the scalar and pseudoscalar contributions cannot be 
compared in the two cases. 

2) The mass splitting can arise from the self-energy diagrams and the tadpole di- 
agrams, but not from the top diagrams. Our conclusion is that the l/A^-corrections 
make the non-singlet heavier than the singlet. 

3) Both in the scalar singlet and non singlet cases the corrections are negative and 
quite sensitive to the value of the UV cut-off. Positive values of the renormalized 
scalar singlet mass are limited to low values of A. We have considered the following 
ranges for the different parameters: 
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250 MeV < M Q < 350 MeV 
0.6 < g A < 0.8 



(21) 

MeV 

500 MeV < A < 2M Q . 



900 MeV < A x < 1100 MeV 



For Mq ~ 250 MeV corrections are big in absolute value (more than 100%). For 
Mq ~ 350 MeV and A ~ 500 MeV we obtain smaller corrections but the mass 
splitting is also smaller. With M Q = 350 MeV, A x = 900 MeV, g A = 0.7 and A = 500 
MeV we find for the singlet M s ~ 500 MeV and for the non-singlet M ns ~ 560 MeV. 



4 Conclusions 

Our work is a first attempt to obtain an order of magnitude as well as the sign of the 
1/N C corrections to the mass of the scalar resonance in the framework of the ENJL 
model. A part from the problem of overlapping divergences, there is a definite way 
of computing these corrections in the full non-bosonized ENJL model. As explained 
in Section 1 we have performed the calculations within the bosonized ENJL version. 

Our main result is the existence at next-to- leading order in 1/N C of a mass- 
splitting between singlet and non-singlet states, making the non-singlet heavier. 

We find that results are extremely sensitive to the values of the UV cut-off A of 
the bosonic loop. Corrections to the singlet mass are negative. However, contrary 
to the mass-splitting, the overall shift of the singlet can be affected by corrections 
coming from higher dimensional fermionic operators. Whereas these corrections are 
about 30% in the vector sector |15|] their importance for the scalar sector is not 
known. 
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Appendix 



Top diagrams 



S 
P 

AO(p°) 
A(V)0(p 2 ) 



i^ii^|2K) 2 r(-2,m 7r ) 



2 /2M|16tt 2 



2r(-2,M A )+4r(-l,M A ) 



. X 3 c A l + 3M 2 
% 2 16tt 2 Mj 

-^3^ W 3^^2Ml (y) r(-2, M A(y) ) 



(22) 



Contributions are the same for singlet and triplet propagator. 

Self-energy diagrams (at q 2 = 0) 
Singlet propagator 



S 
P 

AOijP) x O(p ) 
AO(p°) x 0{p 2 ) 
A(V) 0(p 2 ) x 0(p 2 ) 
A-P 



,Af 1 + 3 
*T 16tt 2 



r(o,M 5 ) 





) 2 1 + 3 [ 


0^ 


' 167T 2 [ 



2m^r(-2, mj - m*r(-l, ra*) + m^r(0, m*) 



^(ca) : 



1 + 3 
16tt 2 



2r(-2, M A ) + r(-l, M A ) + 3r(0, Ma) 



id 



-2ic A J A ^_6Ml 
2 1 + 3 1°Af 4 

167T 2 

2 1 + 3 1 



r(-i,M A )-r(o,M A ) 



2r(-2,M A(lo )-r(-l,M4 ( vo) 



r(o,M 



A(V); 



(23) 



Triplet propagator 

The only possible self-energy diagrams for the triplet propagator contain two differ- 
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ent internal lines, one singlet and one triplet. The contribution is half the contribu- 
tion to the singlet propagator displayed above. 



Tadpole diagrams 
Singlet propagator 



S 
P 



—l Mo 77 

4 s 16tt 2 



r(-i,M s ), 



:4 1) +cf)Ai_y 2m 4 r( _ 2;m7r) 



A(V) O(p ) 
A(V)0(p 2 ) 



„(2) 



1 + 3 



M 2 , 



1{C A(V) + C A(V))J^ IV1 A(V) 

,1 + 3 



2T(-2,M A(v) )+4T(-l,M A{v) : 



i(4lv) + <V))6 T ^2Ml (y) r(-2, M A(y) ) 



(24) 



Notice that in the singlet case there are no contributions from vector vertices of 
order O(p ) like < SV^V» >, because eft + eft = 0. 

Triplet propagator 



5 
P 



o(p° 

A(\/)0(p 2 ) 



-ijt&P - eft + ci 1} + c ft)^- 2 2m 4 n r(-2,m n ) 



) -i(-eft {v) + 3c<g v) + c^V) + c a { v))^ [2r(-2, M A(y) ) + 
- <V) + <V) + 4 1 (V))6 T ^2M4 (y) r(-2, M A(y) ) 



2r(-2,M A(y) )+4r(-l,M A(y) ) 

(25) 
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FIGURE CAPTIONS 



1) Self-energy diagram built by dressing the leading- N c scalar two-point function 

with one loop of linear chains of constituent quark bubbles. 

2) One loop diagrams which give the next-to-leading in 1/N C contribution to the 

scalar propagator: (a) the self-energy diagram, (b) the tadpole diagram and 
(c) the "top" diagram. The solid line is a scalar field, while the dashed line in 
the loop can be any mesonic field. 
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Figure 1: 
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Figure 2: 



15 



